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Indirect interactions of membrane-adsorbed cylinders 

Thomas R. Weikl 

Max-Planck-Institut fiir Kolloid- und Grenzflachenforschung, 14424 Potsdam, Germany 



Abstract. Biological and biomimetic membranes often contain aggregates of embedded or adsorbed macro- 
molecules. In this article, the indirect interactions of cylindrical objects adhering to a planar membrane 
are considered theoretically. The adhesion of the cylinders causes a local perturbation of the equilibrium 
membrane shape, which leads to membrane-mediated interactions. For a planar membrane under lateral 
tension, the interaction is repulsive for a pair of cylinders adhering to the same side of the membrane, and 
attractive for cylinders adhering at opposite membrane sides. For a membrane in an external harmonic 
potential, the interaction of adsorbed cylinders is always attractive and increases if forces perpendicular 
to the membrane act on the cylinders. 



PACS. 87.16.Dg Membranes, bilayers, and vesicles 
and forces, potential energy surfaces for collisions 



34.20-b Interatomic and intermolecular potentials 
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1 Introduction 

Biological membranes consist of a multi-component lipid 
bilayer with a variety of embedded or adsorbed macro- 
molecules such as proteins PE]- In recent years, exper- 
iments revealed a complex lateral architecture of these 
membranes which contain domains or 'rafts' of different 
molecular composition. The domains often serve impor- 
tant biological functions in signaling budding or 
cell adhesion OEI- 

Lateral phase separation and domain formation has 
also been observed for biomimetic membranes, which are 
composed of only a few different molecules. In principal, 
the phase separation may either be caused (i) by a demix- 
ing of the lipid bilayer |7||H], or (ii) by the aggregation 
of embedded or adsorbed macromolecules To un- 

derstand the latter, membrane-mediated interactions be- 
tween macromolecules have been studied intensively. They 
can be divided into static and dynamic interactions. Dy- 
namic, or Casimir-like, interactions arise from the sup- 
pression of membrane shape fluctuations by embedded 
macromolecules such as rigid inclusions |lllll2limill4lll5l 
Iltilll7| or specific receptors or stickers jl8lll9j . Static in- 
teractions are due to a perturbation of the equilibrium 
bilayer structure or equilibrium membrane shape by the 
embedded or adsorbed macromolecules. Examples for such 
macromolecules are trans-membrane proteins ^(3,21,22', 
I^SII^ and adsorbed molecules (2Sj causing a perturbation 
of the equilibrium membrane thickness, as well as con- 
ical or anisotropic inclusions |llll2()ll27ll28ll29lll5IIH0j and 
membrane-anchored polymers |31II32| which cause a per- 
turbation of the equilibrium membrane curvature. 

In this article, the static interactions of parallel cylin- 
drical objects adhering to a planar membrane are con- 
sidered. A single such cylinder has been recently studied 



in Ref. The cylinders are characterized by their ra- 
dius R and a favorable adhesion energy per area U which 
is balanced in equilibrium by the elastic energy of the 
membrane shape deformation. Examples for such objects 
are cylindrical viruses or coated latex cylinders similar to 
the beads considered in Ref. 10 . The membrane shape 
around an adsorbed cylinder is similar to the shape around 
an elongated wedge-shaped inclusion characterized by its 
diameter D and angle a. However, the pair interaction 
energies are different since the contact area of parallel ad- 
sorbed cylinders depends on the distance L between the 
cylinders, while the diameter D and angle a of wedge- 
shaped inclusions are constant. 



2 General model and geometry 

In the absence of adhering objects, the membrane con- 
sidered here is planar and constrained into the x-y plane 
either by a lateral tension a (section 3) or by a harmonic 
potential (section 4). The adhering cylinders are char- 
acterized by their radius R and adhesion energy U per 
area. The cylinders are assumed to be parallel and to be 
much longer than their distance L. The membrane pro- 
file then is approximated by the one-dimensional function 
h{x) measuring the deviation out of the x-y plane. Here, 
X is the cartesian coordinate perpendicular to the cylinder 
axes. The membrane bending energy per area is given by 
^/c/i"(x)^ where n is the bending rigidity and h"{x) is the 
total membrane curvature. 

In order to determine the equilibrium membrane en- 
ergy, we first determine the shape and energy of the mem- 
brane as a function of the contact area with the cylin- 
der(s). For a single cylinder with axis located at x = 0, 
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Fig. 1. Sketch of a membrane with two parallel cylinders ad- 
hering at the same side. 



the contact area is given by |x| < So- For a pair of cylinders 
with axes located at x — ±L, the contact areas are given 
by L — (St < < L + So (see Fig. 1). In equilibrium, the 
system is 'free' to choose its contact area. Therefore, we fi- 
nally minimize the free energy with respect to the contact 
parameters So, S^ and the deviation ho of the cylinders 
out of the x-y plane. In the case of a single cylinder, the 
contact parameter So divides the membrane into the con- 
tact region and the two 'outer' membrane regions with 
> So- In the case of two cylinders, we have two con- 
tact regions, an 'inner' membrane region with |a;| < L — S^ 
between the cylinders, and two 'outer' membrane regions 
with Ixl > L -f So- 



3 Membrane under lateral tension 

In the presence of a lateral tension a, the elastic energy of 
the membrane can be written as 



G 



' (^h"{xf + ^h\xy 



dx 



(1) 



where h{x) is the membrane profile perpendicular to the 
adhering cylinder(s), and k is the bending rigidity. Here, 
i/i'(a;)^ is the local area increase with respect to the x- 
y plane. The membrane profile h{x) outside of the con- 
tact region with the cylinder (s) has to fulfill the Euler- 
Lagrange equation 



h""{x)-eh"{x)=Q 



(2) 



associated with eq. (Q. Here, ^ = \fojK, is a character- 
istic reciprocal length. A general solution of the Euler- 
Lagrange equation l^l is given by 

h{x) = Ci + C2X + C3 exp(-^a::) -|- C4 exp(^a:) (3) 

The tension cr is assumed to be constant here. Strictly 
speaking, this assumption presupposes a membrane area 
reservoir, since the overall area increase of the membrane 
with respect to the reference plane then depends on the 
distance of the cyhnders. 



3.1 Single cylinder 

We first consider a single cylinder adhering to the mem- 
brane in the absence of an external force. The center of the 



cylinder is located at a; = 0, and the membrane adheres 
to the cylinder for —So<x<So- The membrane segment 
in contact with the cylinder has the circular profile 



h{x) ^ ho - \Jb?- -x^ + R~ho 



X 

2R 



(4) 



for So <C R, where R is the cylinder radius. For \x\ > So, 
the profile of the membrane has the form 



h{x) = A + Bexp{-^\x\) 



(5) 



For A = 0, the profile fulfills the boundary condition 
for |a;| ^00. The constants B in eq. ISJ and ho 
in eq. Q can be determined from the conditions of con- 
tinuous profile h{x) and slope h'{x) at a; = So- From the 
latter condition, we obtain B — —So/ {£,R) to first order in 
So- The energy Gz of the adhering membrane with profile 
is the sum of the elastic energy and the adhesion 
energy. To second order in So-) the energy of the ctdliering 
membrane is given by 



- 2So 



U 



(6) 



where U is the adhesion energy per area. The energy of 
the 'outer' membrane segments with |a;| > (5o is 



G„ — 



Si 



i?2 



(7) 



From the equilibrium condition dG/dSo = with G 
Gz + Go, we obtain 



So^ 



K + 2R^U 
2\/aK 



(8) 



The cylinder is bound to the membrane for So > 0, i.e. for 
U < —k/{2R'^). Thus, to obtain a bound state, the ad- 
hesion energy U has to compensate at least the bending 
energy k/{2R^) of the membrane with curvature 1/R at 
the cylinder. Inserting into the total membrane energy 
G = Gz + Go leads to 



G = - 



{k + 2R^Uf 

i?2 



(9) 



The profile of a membrane with a single adsorbed cylinder 
is shown in Fig. [21 




Fig. 2. Minimum energy profile for a single cylinder adhering 
to a membrane under a lateral tension a. The rescaled adhesion 
energy is U = 2UR'^/k, = -2, and the inverse characteristic 



length ^ 
radius. 



/a/K has the value ^ = 1/R where R is the cylinder 
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Fig. 3. Dimensionless interaction energy g{L) — —(1 + 
tanh(^I/)) of two cylinders adhering to the same side of a 
membrane under lateral tension (see eq. USU I. Here, is the 
rescaled distance of the cylinders. 
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Fig. 4. Dimensionless interaction energy g{L) = —(1 + 
coth(^L) of two cylinders adhering to opposite sides of a mem- 
brane under lateral tension (see eq. 121^ '). Here, is the 
rescaled distance of the cylinders. 



3.2 Two cylinders at same membrane side 



3.3 Two cylinders at opposite membrane sides 



We now consider a membrane in contact with two parallel 
cylinders with centers located at a; = ±L. In this section, 
the cylinders adhere at the same side of the membrane. We 
assume first that the membrane is bound to the cylinders 
for L — Si, < \x\ < L -\- 6o, and later determine the contact 
points L — Sl and L + So from a minimization of the energy. 
A general solution for the shape of the 'inner' membrane 
segment with |a;| < L — d,, obeying the symmetry condition 
h{x) — h{—x) is 



h{x) = C + Dcosh{^x) 



(10) 



From the condition of continuous slope h' (x) at x ^ L — S,,, 
we obtain D = —5J{£^Rs'mh{^L)) to first order in S^. To 
second order in S^, the energy of the 'inner' membrane 
segment with — (L — i5t) < a; < i — (5^ is 



15^ 



tanh(CL)i?2 



(11) 



The total energy of the membrane is G = G,, + Go + G^ 
with Go given in eq. and G^ = 2{S, + 6o){U + k/ {2R'^)) 
as in the previous section. The contact parameters and 
So follow from the equilibrium conditions dG/d5i = and 
dG/d5o = 0. We obtain 



5, 



{k + 2R^U) tanh(Ci) 



(12) 



and 5o as in eq. and the interaction energy 



G{L) = 



(K + 2i?2[/)2(i + tanh($L)) 



(13) 



The interaction energy is repulsive, attaining its mini- 
mum G = -{k + 2R^Uf / {2^/^R^) for L^oo. The di- 
mensionless interaction energy g{L) = Ay/aKR'^G{L) / {ti + 
2R^Uf is shown in Fig. El 



In this section, we consider two cylinders adhering at op- 
posite sides of the membrane. We assume that the 'right' 
cylinder with axis at x = L is on top of the membrane, and 
the 'left' cylinder at a: = —L below the membrane. The 
membrane segment with L — 5i^<x<L + 5o adhering to 
the right cylinder then has the profile 



h{x) ^ ho + 



2R 



(14) 



for Si ^ R and 6o ^ R- For symmetry reasons we now 
have h{~x) = —h{x). Therefore, the general solution for 
the profile of the 'inner' membrane segment with |a:| < 
L — (5t is 

h{x) =Cx + Dsinhi^x) (15) 

From the boundary conditions of continuous profile h and 
slope h' at X = L — di,, we obtain 



G 



D = 



^Rho cosh(C£) + S, sinh(^L) 
R{^L coshi^L) ~ sinh(CL)) 

L(5t + Rho 
i?(sinh(CL) - ^Lcosh(^L)) 



(16) 
(17) 



to first order in ho and . To second order in ho and 6^ , the 
energy of the 'inner' membrane segment with < L — 5^ 
is 



G, = 



a[S,R^hl cosh(^L) + 5,{L5, + 2Rho) sinh(gL)] 
i?2[^Lcosh(^L) -sinh(eL)] 



(18) 

As in the previous sections, the energy Go of the two 
'outer' segments is given by eq. (0, and the energy of 
the two adhering membrane segments is G^ — 2{5i + 
So)(U + k/ {2R'^)). Minimizing the total membrane energy 
G = Gz + Gi + Go with respect to S^, So, and ho leads to 



Si = - 



{k + 2R^U) coth(^L) 



2R^U 



2aR 



(19) 
(20) 
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and So as in eq. ((HI, and the attractive interaction energy 

(21) 



G(L) = - 



(K + 2i?2[/)2(i + coth(eL)) 




The dimensionless energy g(L) = A^/anR'^GiL) / {k+2R'^U)^ 
is shown in Fig. 01 



4 Harmonic potential 

In this section, we assume that the membrane is bound 
in an external potential, induced, e.g., by a substrate sup- 
porting the membrane, or an elastic mesh coupled to the 
membrane. In harmonic approximation, the elastic energy 
of the membrane then can be written as 



G 



dx 



(22) 



where m is the harmonic potential strength, and h{x) is 
the membrane profile perpendicular to the cylinders. We 
now neglect the lateral tension a. This is justified as long 
as rj = (rn/ (Ak)Y/^ , the characteristic inverse length for 
a membrane in a harmonic potential, is much larger than 
^ = \/ a / hi. The Ginzburg-Landau equation reads 



h"" (X) + 4.T^^h{x) = 

A general solution of this differential equation is 



h{x) = Ci exp(?7a;) cos(7].x) + C2 exp(ryx) sin(?7a;) 
-I-C3 exp(— 77a;) cos(r7a;) + C4 exp(— 77a;) sin(?7a;) (24) 

4.1 Single cylinder 

The profile of the membrane segment adhering to the 
cylinder for |a;| < 60 is again given by eq. To sec- 
ond order in So and /lo, the energy Gz of this membrane 
segment is the same as in eq. The harmonic potential 
contributes a third-order term to Gz which is neglected 
here. 

For I a; I > Jo, the membrane has the form 

h(x) — Aexp(— Tyjxl) cos(r;|a;|) + _B exp(--?7|a;|) sin(77|x|) 

(25) 

which obeys the boundary condition h{x) —^ for |a;| — *■ 
00. The coefficients A and B follow again from the bound- 
ary conditions of continuous profile h{x) and slope h'{x) at 
x = ±6o- To leading order in do and ho, these coefficients 
are given hy A — ho + 5o/{riR) and B — ho- T^o second 
order in 5o and ho, the energy of the 'outer' membrane 
segments with |a;| > 5o reads 

Go ='^{81 + 27^RSoho + 2TfR^hl) (26) 

The contact point 5o and height ho follow from the equi- 
librium conditions dG/dho ~ and dG/d6o — where 




Fig. 5. Minimum energy profiles for (top) a single adsorbed 
cylinder and (bottom) two cylinders adhering at the same side 
of a membrane in a harmonic potential with strength m. The 
rescaled adhesion energy is U = 2UR^/k = -2, and the inverse 
characteristic length 77 = (m/(4K))^'''* has the value r] = 1/R 
where R is the cylinder radius. 



G = Gz + Go is the total membrane energy. We obtain 

K + 2i?2[/ 



(23) and the energy 



5o = 
ho = 

G=- 



2riK 

■ -t- 2r:^u 



{k + 2R^Uf 
At^kR"^ 



(27) 
(28) 

(29) 



If an external force F is acting on the cylinder per- 
pendicular to the membrane plane, the total energy is 
G — Gz + Go + Fho- Minimizing G now leads to 



FR-2tj{k + 2R^U) 
-FR + ri{K + 2R^U) 

4:1]^KR 



(30) 
(31) 



The cylinder unbinds from the membrane if the contact 
parameter So, and thus the contact area, equals zero. The 
threshold force leading to unbinding is 



Ft = 2r]{K + 2R^U)/R 



4.2 Two cylinders at same membrane side 



(32) 



Let us now consider a membrane with two cylinders ad- 
hering at the same side of the membrane. For symmetry 
reasons, the shape of the 'inner' membrane segment with 
< L — Sc is given by 

h(x) = G sinh(?7a;) sm{hx) + D cosh(?7a;) cos(?7a;) (33) 

The coefficients are again obtained from the boundary 
conditions of continuous profile h{x) and slope h'{x) at 
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section 3.3, we assume that the 'right' cylinder with axis 
at a; = i is above the membrane, and the left cylinder at 
X = —L is below the membrane. The profile of the mem- 
brane has the symmetry h{—x) = —h{x). Therefore, the 
profile of the 'inner' membrane segment with |a;| < L — 
has the general form 

h{x) — C cosh(77a::) sin[/ia;] + D sinh(rya;) cos(?7a;) (41) 

From the conditions of continuous profile h{x) and slope 
h' {x) at a; = ±(L — we obtain to first order in 5^ and 



Fig. 6. Dimensionless interaction energy g{L) — 
— exp{2rjL) / {exp{2r]L) +cos(2ryL) +sin(2r;L)) of two cylinders 
adhering to the same side of a membrane in a harmonic 
potential (see eq. 141)11 *). 



a; = ±(L — (5t). To first order in 6i and ho they read 
2 



C = 



D 



f]R<Pi 
2 



{ [rjRho sm{riL) — cos(77i)] cosh(77L) 

-~r]Rho cos{r]L) sinh(r/L)} (34) 
{ [rjRho cos{rjL) + (5^ sin(ry_L)] sinh(r/_L) 



+r]Rho cosh{i]L) sin(r/L)} 



(35) 



with #1 = sin(27/L) -|-sinh(2?7L). To second order in 5^, and 
hoi the energy of the inner membrane segment is given by 



2rjK 
R^<Pi 



{S^ - 27j^R^hl) cos{2r]L) + {6f + 2rfR^hl) x 



cosh(2?7L) + 2-qR6,ho[ii\nh{2riL) - sva{2riL)] \ (36) 



Minimizing the total membrane energy G — G ^ ~l~ Go ~h G z 
with = 2{5, + 5o){U + k/{2R^)) and Go as given in 
eq. (PSl) with respect to 5^, 5o^ and ho leads to 



<5. = 
5o = 
ho = 



(k + 2R^U){cTq>i2rjL) ~ cos{2r]L)) 
2T]K{exY>{2f]L) + cos{2f]L) + sm{2'qL)) 

{k + 2R'^U){exp{27]L) + cos{2r]L)) 
2riK(exp{2riL) + cos{2riL) + sin{2riL)) 
(k + 2R^U){cyip{2-qL) + cos{2-qL) - s\i\{2-qL)) 



47]2Ki?(exp(2?7L) + cos(2r/i) + sin(2ryL)) 
and the interaction energy 



(37) 
(38) 
(39) 



G(L) 



(k + 2i?2[/)2 exp(2ryL) 



27/Ki?2[exp(2?7L) + cos{2r]L) + sin(2r;L)] 



(40) 

The dimensionless energy 5(L) = 2r]KR'^G{L) / {K+2R'^U f 
is shown in Fig. The global minimum of the interaction 
energy is at Lmin = T^/i^r]) (see also Fig. ISJ. 



4.3 Two cylinders at opposite membrane sides 

In this section, we consider two cylinders adhering at op- 
posite sides of a membrane in a harmonic potential. As in 



G = 



TjR-P: 



■ {riRho[{cosh{riL) cos{rjL) — sinh(77L) sin(77i)] 



-(5t sinh(77L) cos(?7L) } 



(42) 



D = -- 



T]R^2 



{i]Rho[{cosh{T]L) cos{i]L) + sinh(77L) sin(?7L)] 



-S,, cosh{riL) sin(?7L)} 



(43) 



with = sin(277L) — sinh(277L). To second order in 5^ and 
ho, the energy of the inner membrane segment is given by 



G. = 



2r]K 



(6^ - 2ifR^hl) cos(2r;L) - (<5f + 2ri^R^hl) x 



X cosh(2?7L) - 2r]RS,ho [sinh(2?7L) + sin(2r7L)] | (44) 

Minimizing the total membrane energy G — G i ~\~ Gq ~\- G z ; 
with G^ = 2[6,+5o){U+k/[2R^)) and Go given in eq. ^ 
leads to 



5, = - 



(k + 2R^U){exp{2r]L) + cos{2r]L)) 
2riK{exp{2r]L) — cos(2?7_L) — sin(277L)) 
(k + 2i?2f7)(exp(2r/L) - cos(2r/L)) 



2?7K(exp(277L) — cos.{2rjL) — sin(2?7L)) 
^ _ (k -t- 2i?2[/)(exp(2?7L) - cos{2r]L) + sin{2T]L)) 
° ~ 4r72Ki?(cxp(277L) - cos(27?L) - sin(2r7L)) 

and the interaction energy 



(45) 
(46) 
(47) 



G{L) 



{k + 2R^Uf exp(2r/L) 



2?7Ki?2[cxp(27]L) - cos{2riL) - sin(2r;L)] 



(48) 

The dimensionless interaction energy g{L) = 2rjKR^G{L) / {k,+ 
2i?2[/)2 is shown in Fig.[3 

4.4 External force 

In the presence of an external force F acting on the cylin- 
ders perpendicular to membrane plane, the total energy 
for given contact parameters 5^ and So and given ho is 
G = G,(5„ So) + G,(J„ ho) + Go{5o, ho) + 2Fho. Minimiz- 
ing G with respect to So, S^, and ho leads to 



cos{2rjL) 



■ sm 



(2r;L)] + F2[cosh(27/L) + cos(277L)] \ (49) 
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Fig. 7. Dimensionless interaction energy g{L) — 
— exp{2rjL) / {exTp{2TiL) — cos(2ryL) — sin(2r;L)) of two cylinders 
adhering to opposite sides of a membrane in a harmonic 
potential (see eq. I48II ). 



5'min 




Fig. 8. Dimensionless binding energy gmin as a function of the 
rescaled force F for two cylinders adhering at the same side 
of a membrane in a harmonic potential. The absolute value of 
the binding energy ^niin increases under an applied force F, 
irrespective of the sign of F. 



with ^3 = exp(277L)-|-cos(2?7L)+sin(277L) and the rescaled 
force F = -FR/[r]{K + 2R^U)]. 

Fig. |S1 shows the dimensionless binding energy of the 
cylinders 5„,in = (G(L„,in)-G(oo)) x (27?Ki?2)/(^+2i?2[/)2 

as a function of the rescaled force F. Here, Lmin is the 
cylinder separation at which the interaction energy G{L) 
has its minimum. The binding energy of the two cylinders 
increases irrespective of the sign of the force, i.e. both 
for perpendicular forces pulling on the cylinders, and for 
forces pushing the cylinders into the membrane. 



5 Discussion 

We have considered the indirect interactions of parallel 
cylindrical objects adhering to a planar membrane. The 
interactions arise from the perturbations of the equilib- 
rium membrane shape caused by the cylinders. The cylin- 
ders are assumed to be much longer than their distance 
L, and the membrane profile is approximated by a one- 
dimensional function h(x). Here, a; is a cartesian coordi- 
nate perpendicular to the cylinder axes. 

For a pair of cylinders adhering to a membrane under 
lateral tension, the indirect interaction is repulsive if the 



cylinders adhere to the same membrane side, and attrac- 
tive for two cylinders adhering at opposite membrane sides 
(see section 3). For a membrane in a harmonic potential, 
the interaction of the cylinders is always attractive (see 
section 4). The latter situation applies, in first approxi- 
mation, to supported or bound membranes. 

These results are intuitively plausible: The curvature 
h" {x) around a single cylinder adhering to the top of a 
membrane under lateral tension is negative (see eq. (O and 
Fig. [21 . This means that the membrane is 'curved away' 
from a second cylinder approaching the membrane from 
the same side, but 'curved towards' a cylinder approach- 
ing the membrane from the opposite side. The negative 
curvature makes it more difficult for a cylinder to adhere 
at the same side close to the first cylinder, but facilitates 
the adhesion of a second cylinder at the opposite side. In 
contrast, the profile around a single cylinder adhering to 
a membrane in a harmonic potential is an exponentially 
damped, oscillating profile, with both negatively and pos- 
itively curved segments (see eq. (I25|) and Fig.jSJ. 

For a membrane under lateral tension, the profile around 
a single adhering cylinder with axis at a; = is propor- 
tional to exp(— ^|a;|) (see section 3.1). Here, ^ = \/ a j k is 
the characteristic reciprocal length of a membrane with 
tension a. The same profile is obtained around a sin- 
gle wedge-shaped inclusion in a tense membrane (see ap- 
pendix B.l). Similar membrane profiles are also obtained 
for a membrane with two cylinders or two wedge-shaped 
inclusions. However, the contact region of the adsorbed 
cylinders depends on the distance L of the cylinders, whereas 
the angle a and width D of the wedge-shaped inclusions 
are constant. Therefore, the interaction energy as a func- 
tion of the distance L is different for adsorbed cylinders 
and for wedge-shaped inclusions (see eqs. ((T^ . (|^ . ifTUi) . 
and JHOI))- For the same reason the interaction energy of 
adhering spheres should be different from the interaction 
energy of conical inclusions. 

The strength of the cylinder interactions can be esti- 
mated from the energetic prefactors. For a membrane un- 
der lateral tension, the prefactor (k -|- 2i?^?7)^/ (J^^faK,F(^) 
of the interaction energies (|13|l and (|21|l can be written as 
{5o/ RYy/o^, where 5o is the contact point at large separa- 
tions L (see eq. inj- For {So/R) < 0.5, the typical bending 
rigidity k = 10~^^ J and membrane tensions up to 10"'^ 
N/m, the prefactor can attain values of up to 0.5 fcsT per 
nanometer length of the cylinders, leading to large inter- 
action energies for colloidal cylinders with lengths of sev- 
eral hundred nanometers. For a membrane in an external 
potential, the prefactor {k + 2R^UY / {2riKR^) of the inter- 
action energies (|40() and H48() can be written in the form 
2{5o/Rfvn with 5o given by eq. For 6o/R = 0.5, a 

typical length scale 1/r] = 50 nm of about 10 times the 
membrane thickness, and n = 10~^^ J, this prefactor has 
the value 0.25 fc^T/nm. The binding energy for two cylin- 
ders adhering at the same side then is about -0.025 fcsT 
per nanometer length of the cylinders (see Fig.l^j), leading, 
e.g., to a binding energy of 10 fc^T for colloidal cylinders 
with a length of 400 nm. These estimates for the static 
interactions of the colloidal membrane-adsorbed cylinders 
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are clearly larger than the loss of orientational free energy 
of order fc^T per cylinder, which is caused by the parallel 
alignment. The interactions therefore can lead to stable 
pairs or bundles of such cylinders. However, for a mem- 
brane dominated by lateral tension, these conformations 
are only stable if the cylinders are adsorbed alternately on 
opposite membrane sides. 

Attractive interactions and aggregation of spherical 
particles adsorbed on vesicles has been experimentally ob- 
served by Koltover, Radler, and Safinya PHI- The parti- 
cles do not aggregate in solution. Therefore, the attractive 
interactions are probably indirect, i.e. mediated by the 
membrane. Since the particles have a radius of 0.9 fim, the 
interactions are most likely due to the rather long-ranged 
fluctuation- induced interactions, or interactions from per- 
turbations of the equilibrium membrane shape as consid- 
ered here. Interactions due to membrane thickness pertur- 
bations typically decay over a few nanometers, a length- 
scale comparable to the membrane thickness 20,21,22] 
I23II24II25| . The fluctuation-induced pair interactions are 
always attractive, but rather weak compared to thermal 
energies |lllll2lll8lll4lll5lll6j . Depending on the induced 
deformation, the static interactions considered in this pa- 
per can be strong, see above. Exploring the latter as a pos- 
sible explanation for the observations made by Koltover 
et al. requires an extension of the calculations presented 
here to beads and to nonplanar membranes or vesicles. 



the elastic energy of the two outer membrane segments is 



ClL + Clai - e-2eA) ^ clae^^"" - 1) 
+2C2C3{e-^^ - 1) + 2C2Ci{e^^ ~ 1) 



(53) 



For given ho and 5o, the analytical expressions for the four 
parameters Ci to C4 are rather lengthy and not shown 
here. However, if we minimize the total free energy G = 
Gz + Go with Gz as in eq. © with respect to ho and So, 
we get 



C2 
C3 

Ga 



{k + 2R^U) 
~ 2aRcosh{£,A) 


(K + 2i?2[/)(tanh(^yl) + l) 
(K + 2i?2[/)(tanh(^/l) - 1) 



(54) 
(55) 
(56) 

(57) 



For yl ^ 00, we obtain Ci -> 0, C3 -> {k + 2R^U) / {2aR), 
and C4 exp(^a;) — > 0. In the limit of infinite membrane 
size, we thus recover the equilibrium shape h{x) = [n + 
2R'^U)exp{-£,x)/{2aR) from section 3.1. 
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A Single cylinder adhering to a finite 
membrane under lateral tension 

In this appendix, we consider a single cylinder adhering to 
a finite membrane under lateral tension. In the limit of in- 
finite membrane size, we will obtain and thus corroborate 
the profile and energy derived in section 3.1. 

As in section 3.1, the membrane segment in contact 
with the cylinder for —do < x < 60 has the circular profile 
(0J). The profile h{x) of the nonadhering finite membrane 
segment for < a; < yl has the general form (O for x > 0. 
For symmetry reasons, we have h{—x) = h{x). The four 
parameters Ci to C4 in © can be determined from the 
four boundary conditions: 

HSo)^ho + $-, hiA)=0 (50) 

h'iSo)^^-^, h'{A)=0 (51) 

(52) 

Here, we again assume So ^ R, and determine Ci to C4 
to first order in So and ho- The second order in So and ho. 



B Wedge-shaped inclusions 

B.l Single inclusion 

A wedge-shaped inclusion is characterized by its diameter 
D and angle a. If the center of the inclusion is located at 
a; = 0, a single inclusion imposes the boundary condition 

h'{±D/2) = ±a (58) 

For a membrane under lateral tension, the profile around 
a single inclusion is given by 

h{x) = A + B cxp(-^| xl) (59) 

as in the case of a single adhering cylinder (see section 
2.1). From the boundary condition (|58|l . one obtains 

S = -aexp(^L>/2)/e (60) 

The coefficient A in H59|) is arbitrary and defines the ref- 
erence plane. The elastic energy of the membrane around 
the inclusion is given by 

G = a^VoK (61) 

B.2 Two inclusions with equal orientation 

We now consider two inclusions of diameter D with 'inner' 
edges located at a; = ±i, and 'outer' edges located at 
X = zk{L + D). In this section, the two inclusions have the 
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same orientation. At their edges, the inclusions impose the 
bomidary conditions 



h'{±L) = ±{P-a) 
h'{±{L + D)) = ±{f3 + a) 



(62) 
(63) 



on the membrane where (3 is the tilt angle of the inclusions. 
The membrane profile is symmetric with h{—x) = h{x), 
and decays to zero for 2: oo. The profile of the 'inner' 
membrane segment between the inclusions is given by 

h{x) ^ C + D coshi^x) (64) 

From the boundary condition (|62l) . one obtains 

P — a 



D 



esinh(ei) 



(65) 



The energy of the 'inner' membrane segment then is given 

by 

G, = (/3 - a)2 Va^coth(ei) (66) 

The 'outer' membrane segments for \x\ > (L + D) have 
again the shape (|60|l . The boundary condition (|63|l is ful- 
filled for 

B^~if3 + a)eMaL + D))/C (67) 

and the elastic energy of the two 'outer segments is 

Go = {a + fJ)^^ (68) 

As in the previous section, the coefficient A in H59|) is arbi- 
trary. Minimizing the total membrane energy G ~ G^ + Go 
with respect to the tilt angle (3 of the inclusions leads to 
the equilibrium value 



13 = aexp(-2Ci) 
and the interaction energy 



G{L) = 2a^V^(l + exp(-2^L)) 



(69) 



(70) 



B.3 Two inclusions with opposite orientation 

If the two inclusions have opposite orientation with respect 
to the membrane plane, the boundary conditions at the 
inclusion edges are 



and 



h'{±L) ^{I3~a) 
h'{±{L + D)) = {P + a) 



h{±L) = ±ho 



(71) 
(72) 



(73) 



where ho is the deviation of the inclusions out of the ref- 
erence plane, and (3 is again the tilt angle. We now have 
the symmetry h{—x) = —h{x). The profile of the 'inner' 
membrane segment between the inclusions is given by 



From the boundary conditions above, we obtain 
^ho cosh(^L) + {a- f3) sinh(^L) 



C = 



D = — 



cosh(^L) - sinh(^L) 
ho + {a~ f3)L 



cosh(^L) - sinh(^L) 
The energy of the 'inner' membrane segment is 



(75) 
(76) 



G, = - {2C^L + 4:GD sinh(CL) + D^^ sinh(2^L)) (77) 

with C and D as given above. As in the previous section, 
the two outer membrane segments have the energy (|68ll . 

Minimizing the total membrane energy G — G^ + Go 
with respect to the tilt angle (3 and deviation of the 
inclusions leads to the equilibrium values 



(3 = ~aexp{-2^L) 
ho - -a{l - exp(-2ei))/f 

and the attractive interaction energy 



G{L) = 2a2v^(l - exp(-2^L)) 



(78) 
(79) 



(80) 



h{x) ^Gx + D sinh(^a;) 



(74) 
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